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We determine the neutrino spectra arising from low-mass (4 — 10 GeV) dark matter annihilating 
in the sun. We also determine the low-mass dark matter capture rates (element by element in the 
sun), assuming dark matter interacts either through elastic contact interactions, elastic long-range 
interactions, or inelastic contact interactions. These are the non-detector-specific data needed for 
determining the sensitivity of a neutrino detector to dark matter annihilating in the sun. As an 
application, we estimate the sensitivity of a one kiloton liquid scintillation neutrino detector (such 
, ^ as KamLAND) and LBNE (LAr-based) to low-mass dark matter with long-range interactions and 

^ compare this to the expected CDMS sensitivity. It is found that KamLAND 's sensitivity can 

^ exceed that obtainable from the current CDMS data set by up to two orders of magnitude. 
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I. INTRODUCTION 



There has been great recent interest in low-mass dark matter {mx ~ (9(10) GeV) as a pos- 
sible explanation for the event rates observed at the DAMA CoGeNT [2] and CRESST [5] 
experiments. However, negative results from XENONlO/100 [U E] and CDMS [BHS] have 
created some tension with these positive signals. There has been much discussion in the 
literature regarding experimental issues with all of this data, but the explanation is far 
from clear P-ITT]. It is important to find new tests of this data, particularly tests that do 
not suffer from the same difficulties as direct detection experiments at low recoil energy. 
Collider, gamma-ray, and anti-proton flux search strategies have been employed for this 
purpose [IHHIT]. 

A promising avenue for cross-checking the low-mass direct detection data is with neutrino 
detectors , which search for the flux of neutrinos arising from dark matter annihilation 

in the core of the sun. Dark matter is captured by the sun via scattering off solar nuclei. 
If the sun is in equilibrium, then the dark matter capture rate determines the dark matter 
annihilation rate. The neutrino flux thus constrains the dark matter-nucleus scattering cross- 
section, and allows neutrino detectors to cross-check direct detection experiments without 
some of the particle physics and astrophysics uncertainties that plague other types of indirect 
detection searches. Moreover, the 0{GeV) neutrinos produced from the annihilation of 
low-mass dark matter are easily distinguishable at water Cherenkov, liquid argon or liquid 
scintillator-based neutrino detectors. 

In order to determine the event rate expected at a neutrino detector, one must calculate 
the rate at which dark matter is captured by the sun due to scattering and the neutrino 
spectrum arising from the decay of the dark matter annihilation products. Numerical pack- 
ages such as DarkSUSY [35] are commonly used for obtaining these rates, which are deter- 
mined from numerical simulations. However, the required simulations have not been run 
for masses less than 10 GeV. Moreover, the capture rate has only been calculated assuming 
that dark matter scatters elastically via a contact interaction. Recent models for reconciling 
the low-mass direct detection data have included the possibility of dark matter scattering 
via long-range forces ED E2] , inelastic scattering jSHHll] , and dark matter interactions 
which are isospin- violating [501 ESI SSHIZ] • To determine the sensitivity of neutrino detectors 
to these models, new capture rate calculations must be performed. 

In this work, we calculate the required neutrino spectra and capture rates for low-mass 
dark matter in the sun. We not only consider the capture rate for elastic contact scattering, 
but also for inelastic dark matter and for models in which the dark matter-nucleon interaction 
is mediated by a low-mass particle. We also determine the regions of parameter-space of 
these models for which dark matter in the sun is in equilibrium. 

As an application of these techniques, we consider the sensitivity of a 1 kT liquid scintilla- 
tion detector with 2135 live days of data (roughly the same exposure as KamLAND) to dark 
matter with long-range interactions. We find that the sensitivity of neutrino detectors to 
dark matter with long-range interactions is enhanced, because typical scatters off low-mass 
targets (such as the hydrogen and helium in the sun) involve small momentum transfers, 
yielding enhanced scattering cross-sections in models where the mediating particles are light. 
This implies the existence of an entire class of dark matter models for which current neutrino 
experiments can provide the leading sensitivity. 

In section II, we review the general formalism for dark matter searches using neutrino 
detectors. In section HI, we describe the details of the computation of the neutrino spec- 
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trum arising from the annihilation of low-mass dark matter. In section IV, we describe the 
calculation of the dark matter capture rate in the sun, assuming dark matter with either 
elastic contact interactions, elastic long-range interactions, or inelastic contact interactions. 
In section V, we describe the range of circumstances under which low-mass dark matter is 
in equilibrium in the sun. In section VI, as an example, we apply these techniques to deter- 
mine the sensitivity of a 1 kT LS neutrino detector to low-mass dark matter with long-range 
interactions. We conclude in section VII. 



II. OVERVIEW OF DARK MATTER DETECTION VIA NEUTRINOS 

The rate of charged lepton events at a neutrino detector can be written as 

where d ~ 1 AU is the earth-sun distance, z = E^/rrix, f] is the nucleon number density of 
the detector (including the earth around the detector, in the case of a search for throughgoing 
muons), and e is the efficiency for a neutrino charged-current interaction to produce a charged 
lepton which will pass the detector analysis cuts. Bf is the branching fraction to each 
dark matter annihilation product /, and dNj^^Jdz is the differential neutrino spectrum per 
annihilation to each final state, for each (anti-)neutrino fiavor. Ta is the total dark matter 
annihilation rate, and (Tu^-n is the (anti-)neutrino-nucleon scattering cross-section. For dark 
matter in the 4 — 10 GeV range, most of the charged leptons produced in a reasonably-sized 
detector will be fully-contained, with the vertex where the lepton is produced and the end 
of the track both within the fiducial volume of the detector. As a result, we will focus on 
fully-contained lepton events. 

The quantities r] and e are specific to the geometry and construction of the detector. For 
1 GeV < Ey <1 TeV, the a^^-N can be approximated as [48] 

(T„_„ ~ 4.51 X 10"^ X z {ttix/ GeV) pb, 

(T^_„, ~ 8.81 X 10"^ X z{mx/ GeV) pb, 

ap_p ~ 3.99 X 10"^ X z (m^/ GeV) pb, 

ap_„ ~ 2.50 X 10"^ X z {nix/ GeV) pb. (2) 

The two remaining quantities we will need to compute are and dNf/dz. 
For a search for fully-contained charged leptons, we can then write 

«=|5^x/'''^E^/^'-x^.//.W. (3) 

where 

AeffXz) = (Ty^^N{Ey = mx) j dVri{r) X e{r,z). (4) 

All detector-specific information is encoded in the effective area, Aef/X^)- 

If ^ 4 GeV, the effect of dark matter evaporation can be important [M]. In 
that case, dark matter in the sun's core is significantly depleted by evaporation, and the 
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total annihilation rate is relatively small, implying that constraints on dark matter from the 
neutrino flux will be weak. We focus on the regime mx > 4 GeV, for which dark matter in 
the sun can only be depleted through annihilation. If the sun is in equilibrium, we then find 
that is related to the dark matter capture rate, Tc, by the relation Tc = 2Ta- Since 
is determined by the dark matter-nucleus scattering cross-section, the above relation allows 
one to translate neutrino flux bounds into bounds on the dark matter-nucleon scattering 
cross-section. 

III. NEUTRINO SPECTRA 

For low- mass dark matter {mx < 20 GeV), the only dark matter annihilation products 
relevant for neutrino searches are heavy quarks (6, c), r, i'e,ti,T and the gluon (g). Muons 
and light quarks {u, d and s) will tend to stop within the sun before they decay |19]. The 
resulting neutrinos are thus very soft (though they also can potentially be used for for dark 
matter searches [50]). 

We have calculated the differential neutrino spectrum per annihilation for the rele- 
vant channels using the DarkSUSY/WimpSim/NuSigma/Pythia [351 EH 152] package. Nu- 
merical simulations were run on the Hawaii Open Supercomputing Center (HOSC) com- 
puting cluster. 10^ annihilations were simulated for each annihilation channel and dark 
matter mass in the range 4 — 10 GeV (in increments of 2 GeV). Representative spec- 
tra are plotted in Appendix A, and all of the original data files are available at |http : | 
77wwwTphysThaijaii.edu/~superk/post/spectrum. We present the neutrino spectrum at 
a distance 1 AU from the sun, including the effects of hadronization and decay of the an- 
nihilation products at injection, matter effects as the neutrinos propagate through the sun 
(including tau-regeneration) and vacuum oscillations. The neutrino oscillation parameters 
were chosen to be 

= 33.2°, 

023 = 45°, 
Am^i = 8.1 X 10"^ eV^ 
Am^2 = 2.2 X 10"^ eV^ 

^13 = 10°, (5) 

assuming a normal hierarchy. This choice is consistent with recent exciting data from the 
Daya Bay experiment [53] indicating ^13 ~ 9°. Data files for the choice ^13 = 0° are also 
available online. For this choice, the change in the neutrino spectrum is relatively small 
unless dark matter annihilates directly to neutrinos in a flavor-dependent way. 

In the case where dark matter annihilates to 6-quarks, annihilation can only proceed if 
the dark matter mass is larger than the mass of the 6-hadron which is produced. Moreover, 
a 6-quark will lose ~ 27% of its energy during hadronization ^\ IM] . Thus, the neutrino 
spectrum arising from annihilation to 6-quarks is simulated only for mx > 6 GeV. In the 
case where dark matter annihilates to rf, the neutrino spectrum has been computed by 
averaging over helicities. Some dark matter candidates will preferentially decay to certain 
helicities, which can have a significant effect on the injected neutrino spectrum |[55]| . 

The neutrino spectrum at 1 AU is equivalent to the spectrum of downward-going neutrinos 
at the detector, averaged over the year. This spectrum determines the rate of downward 
going fully-contained charged leptons. In the case where the charged lepton is upward going 
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through the earth, one should also include oscillation and matter effects as the neutrinos 
pass through the earth. These effects depend on the location of the detector and can be 
determined by inputting the neutrino spectrum at 1 AU into the "WimpEvent" program, 
with the location of the detector specified. 

IV. CAPTURE RATE 

The capture rate can be computed following the analysis of and we follow that 
notation. A dark matter particle in the halo has velocity u, given by a distribution f{u) 
obeying J duf{u) = rjx- Here, rjx is the dark matter number density in the halo. When 
a dark matter particle is at distance r from the core of the sun, it will have velocity w = 
+ f^, where v{r) is the escape velocity from the sun at radius r. Thus, a dark matter- 
nucleus scatter will result in dark matter capture if the dark matter scatters from velocity 
w to a velocity < v. 

More generally, however, 3-body interactions can drastically affect the capture rate. As 
a conservative estimate, one can choose to count as "captured" only dark matter particles 
which are kinematically constrained to orbits with maximum radius tq, often taken to be 
the radius of Jupiter's orbit. The velocity to escape from position r < within the sun to 
radius is denoted by fe(r), and is given by the relation fe(r)^ = v{r)^ — f (ro)^. 

The rate for dark matter to be captured in any differential solar volume by scattering off 
an element with atomic number Z can then be written as 



where Umin,max are the minimum/maximum dark matter velocities in the halo such that 
a nuclear scatter at position r resulting in dark matter capture is kinematically possible. 

(w) is the rate per unit time at which a dark matter particle with velocity w will scatter 
to velocity < Ve{r), and is given by the expression 



where tjq is the number density of the sun, Emin is the minimum recoil energy needed for 
capture, and Emax is the maximum recoil energy that is kinematically allowed. da^'^/dEn 
is the differential cross-section for dark matter to scatter off a nucleus with Z protons and 
A nucleons. 

For any dark matter model, the quantities that must be known to compute the cap- 
ture rate are da^'^^/dE^, Emin,max, and Umin,max- Given these quantities, we can compute 
the capture rate numerically using the DarkSUSY code (appropriately modified), with its 
standard assumptions about solar composition. 




(6) 




(7) 
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A. Elastic Contact Interactions 



The most commonly used assumption is that dark matter interacts with nuclei via elastic, 
isospin-invariant, contact interactions. In this case, we have 
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where fi = mx/mA and /x± = (/i ± l)/2. are the relative strengths of dark matter 
coupling to protons and neutrons, respectively, is the dark matter-proton scattering 
cross-section, and Fa{Er) is the nuclear form factor. To match the assumptions used in 
DarkSUSY, we will assume a Gaussian form factor \Fa{Er)\'^ = exp[—Eji/Eo\, where Eq = 
3h^ /2mARA- Ra = [0.3 + 0.91(771^/ GeV)^/^] fm is taken as the nuclear radius. 

Following default assumptions in DarkSUSY [351 ES], we have assumed a Maxwell- 
Boltzmann velocity distribution for dark matter in the Galactic halo of the form 



fhalo{u) = rjx 



u 



(9) 



where v is the three-dimensional velocity dispersion, which we have set to v = 270 km/s. 
This velocity distribution is truncated at the galactic escape velocity, Vesc- The velocity 
distribution seen by an observer moving through the halo with velocity is then 
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The step function imposes the condition that f{u) = for m > t>=K + Vesc- If we ignore the 
truncation at the galactic escape velocity, this reduces to the expression 



f{u) = fhaio{u)e 



-3t)2/2s2_^ 



sinh 



3uw* 



We take = 220 km/s to be the velocity of the sun through the halo. 
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FIG. 1: Solar dark matter capture rates for various elements in the sun, assuming isospin-invariant 
elastic contact interactions with a'^ = 10~^ pb. 



For isospin-invariant interactions, one assumes fn/fp = 1- For isospin- violating dark 
matter, the capture rate from scattering by each element is scaled by a factor fp) {A— 

Z)Y/A^. To facilitate this rescaling in the case of generic isospin-violating interactions, we 
plot the capture rate for each of the main elements in the sun separately. 

We have plotted in fig. [T] the capture rates for elastic contact interactions if one requires 
dark matter to be captured to within the radius of Jupiter's orbit. If the presence of Jupiter 
is neglected, capture rates change by less than 1% in the case of elastic contact interactions. 



B. Elastic Long- Range Interactions 

If dark matter interacts with nuclei via long-range interactions (or, equivalently, via a t- 
channel interaction with a mediator with a mass much smaller than the momentum transfer), 
then the differential scattering cross-section will have a different form. The quantum matrix 
element scales as oc l/(g^ — M^) where is the mass of the mediating particle and q is 
the momentum transfer; for <^ q"^, the differential cross-section scales as q~^. For such 
a model, it would not make sense to parameterize the differential scattering cross-section in 
terms of o"p since, as with Rutherford scattering, the total cross-section is infinite. We may 
instead write 
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FIG. 2: Solar dark matter capture rates for various elements in the sun, assuming isospin-invariant 
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elastic long-range interactions (as in eq. 13 ) with C 
be confined within the orbit of Jupiter. 



10 . Captured dark matter is required to 



where /Xp = rnxfrip/ {mx + Tn^) is the dark matter-proton reduced mass. C is a constant 
which defines the size of the differential scattering cross-section in terms of the proton 
charge; if gx,p are the strengths with which the mediator couples to the dark matter and a 
proton, respectively, then C = gj^Qp/e^. Using the Gaussian form factor defined above, the 
integrated differential scattering cross-section takes the simple form 
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(14) 



The kinematics of the scattering process are the same as in the case of an elastic contact 
interaction, and thus Emin,max and Umin,max are the same as in equation (|8|. The capture 
rates for dark matter with long-range interactions are plotted in fig. |2] (for mx = 4 — 10 GeV) 
and fig. [s] (for mx = 10 — 1000 GeV), again assuming that captured dark matter must be 
confined to an orbit inside Jupiter's. 

One should note that, in the case of long-range interactions, it is necessary to assume 
that captured dark matter be confined to an orbit within some finite radius tq; without this 
assumption, the capture rate would be infinite. The origin of this divergence is easily un- 
derstood to arise from the low-velocity tail of the Maxwell-Boltzmann velocity distribution. 
Near u = 0, dark matter far from the sun has a very small kinetic energy. As a result, even 
scattering interactions yielding very small recoil energies can result in a dark matter particle 
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FIG. 3: Solar dark matter capture rates for various elements in the sun, assuming isospin-invariant 
elastic long-range interactions (as in eq. 13) with C = 10"^". Captured dark matter is required to be 
confined within the orbit of Jupiter. Plotted for dark matter with mass in the range 10 — 1000 GeV. 



being captured, (i. e., having negative total energy). Since the differential scattering cross- 
section diverges at small recoil energy, the total capture rate diverges. This simply reflects 
the fact that it is not physically sensible to think of dark matter as captured if confined 
to an orbit of very large radius. It is most sensible to count as captured only dark matter 
confined to orbits that lie within Jupiter's orbit. 

Note that we are not including the possibility of capture due to multiple scattering. For 
many models, these effects can significantly enhance the dark matter capture rate, especially 
in the case of long-range interactions. The dark matter capture rate may be much larger if 
dark matter scattering exhibits Sommerfeld enhancement. But this depends on the details 
of the model, including the nature of dark matter interactions with electrons and possible 
3-body effects. These issues may be relevant for specific models but are beyond the scope 
of this work. 



C. Inelastic Contact Interactions 

One may also consider the case where dark matter scatters inelastically off nuclei, via 
the process XA — )■ X'A. We will consider the case with 6mx = ^x' — > 0. In this 
case, the scattering matrix element will only change by subleading 0{5mx/'fnx) terms, but 
the kinematics of the scattering process can change dramatically. It is easiest to consider 
this process in the center-of-mass frame. We then find pj — p'j ~ 2mrS'mx, where rrir = 
^xf^A/i^x + fnA) is the reduced mass, and Pi = rrirW and pf are the spatial momenta 
of the incoming X and outgoing X', respectively, in the center-of-mass frame. The phase 
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space factor of the differential scattering cross-section is directly proportional to the outgoing 
momenta. 

Again, it is not appropriate to express the dark matter-nucleus inelastic scattering cross- 
section in terms of the dark matter-proton scattering cross-section, since there exist kine- 
matic regions where dark matter-proton inelastic scattering is impossible, though dark mat- 
ter can scatter off other nuclei. Using the fact that the recoil energy can be written as 
= {2mA)~^{Pi + p} — 2piPf cosOcm), we can write 
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where m\m\[Z + {A — x / is the squared dark matter-nucleus matrix 

element (summed over final spins and averaged over initial spins). / is roughly constant for 
different elements (up to 0{6mx/mx) corrections), so it makes sense to present bounds on 
inelastic dark matter in terms of this quantity. We also find 
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where the first term is the minimum recoil energy kinematically possible in two-body inelastic 
scattering, and the second term is the minimum recoil energy in a process where the outgoing 
dark matter particle is slower than Vg, the velocity to escape to radius tq. 
Umax is determined by the constraint Emax > Emm, yielding 

" 26mx fi^ (1 + yu')t^(ro)' 

7/ = —7; — 
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Finally, we have 
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(18) 



(19) 



because inelastic scattering is only kinematically possible for 5mx < (l/2)m.rW^. This 
implies that, for mx ^ 10 GeV one need only consider models with 5mx ^ (9(10 — 100) keV. 

The capture rates for inelastic dark matter with contact interactions and 5mx = 
10, 30, 50 keV are plotted in fig|4j It is interesting to note that, as 5mx increases, the rate 
of capture arising from scattering off light elements vanishes. This is because m^ is smallest 
for light elements, implying that they have the smallest maximum value of 5mx such that 
inelastic scattering is kinematically allowed. It is also worth noting that inelastic scattering 
for low-mass dark matter is kinematically allowed in the sun for larger 5mx than in the 
Earth, because dark matter within the sun has gained kinetic energy from gravitational in- 
fall. As a result, even low-mass inelastic dark matter with 5mx ~ 50 keV can be potentially 
probed by neutrino detectors. 
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FIG. 4: Solar dark matter capture rates for various elements in the sun, assuming isospin-invariant 
inelastic elastic contact interactions with //327r = 10~^ pb GeV~^. The three panels are for 
5mx = 10, 30, and 50 keV, as labelled. 



V. EQUILIBRIUM 

If the effects of WIMP evaporation are negligible, the equilibration time Tq for dark 
matter in the sun can be written as [Ml l57] 

where ^0 ~ 4.5 x 10^ yr is the age of the solar system. 

For the case of elastic contact interactions (assumed to be isospin-invariant, and either 
spin- independent or spin-dependent), we plot in the top panel of fig. [s] the minimum 
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FIG. 5: Minimum (top panel) or C (bottom panel) required for dark matter to be in equilibrium 
in the sun, assuming either elastic contact (spin-dependent or spin-independent) or long-range 
interactions, respectively. We assume that the dark matter annihilation cross-section is given by 
(av) = 1 pb. 

required for the sun to currently be in equilibrium, assuming that the total dark matter 
annihilation cross-section is given by (erf) = 1 pb. Note that, for the case of IVDM with 
spin-independent interactions, the required for the sun to be in equilibrium would lie 
between that required for spin-dependent scattering and that required for isospin-invariant 
spin- independent scattering. Note that the equilibration time scales as {a^{av))~^^'^. IVDM 
with spin-independent interactions and fn/ fp ~ —0.7 could be consistent with the data of 
DAMA, CoGeNT and XENONlO/100 if (rgj ~ lO^^ pb [IHl HE]. Such dark matter can be 
in equilibrium in the sun even if (av) ~ 10~^ pb. This implies that, if the IVDM candidate 
is a thermal relic, then it can currently be in equilibrium in the sun even if almost all of 
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the annihilation cross-section at freeze-out was due to p-wave interactions (suppressed at 
current times), with only a negligible amount due to s-wave interactions. 

Similarly, for the case of long-range interactions, we plot in the bottom panel of fig. [5] the 
minimum C required for the sun to be in equilibrium (assuming (cxf) = 1 pb). 

VI. AN APPLICATION TO NEUTRINO DETECTORS 

We now consider an application of these tools to a specific detector. We will focus 
on the case of long-range interactions, because neutrino detectors are expected to have a 
major advantage over direct detection experiments in this instance. For both direct de- 
tection experiments and neutrino searches, the measured event rate will be proportional 
to the dark matter-nucleus scattering cross-section. For the case of long-range interac- 
tions, the integrated dark matter-nucleus scattering cross-section is roughly proportional to 
Cm'^^ E^lj^E~l^^, where uia is the mass of the nucleus and Emin is the minimum nuclear 
recoil energy which one can measure. For a direct detection experiment, Emin is the recoil 
energy threshold of the experiment, and is typically of order 2 — 10 keV. For germanium- 
based experiments (such as CDMS and CoGeNT), ~ 72 nip, while for xenon-based 
experiments ~ 130 rrip. For neutrino searches, -Emm is the minimum recoil energy 
such that dark matter is captured and can annihilate in the core of the sun. We thus find 
Emin ~ (l/2)mx'U^ ~ 2 — 5 keV, for dark matter in the mass range considered here. How- 
ever, mA is the mass of the target nucleus in the sun and is very small for some elements 
that contribute significantly to capture in the sun. For example, hydrogen contributes ~ 3% 
of the dark matter capture rate for low-mass dark matter, and tjih = mp. So one can expect 
the sensitivity of a neutrino search for low-mass dark matter with long range interactions to 
be significantly enhanced (~ 10^ — 10^) compared to direct detection experiments. 

We compare the sensitivity of liquid scintillation (LS) neutrino detectors to that of CDMS. 
It was shown in [SS] that liquid scintillation neutrino-detectors can determine the flavor and 
direction of leptons produced by a charged-current interaction using the timing of the first 
photons which reach the photomultiplier tubes. We will focus on a search for electron 
neutrinos producing fully-contained electron/positron events. An advantage of this strategy 
is that the atmospheric electron neutrino background is significantly smaller than that of mu 
neutrinos. It was estimated that liquid scintillation neutrino detectors can provide almost 
absolute lepton flavor discrimination, and electrons of the energy range we consider can be 
measured with an angular resolution < 1°. It was also estimated that the neutrino energy 
could be determined (from the energy and direction of the produced charged lepton, as well 
as total energy deposition) with a resolution ~ 1 — 3%. 

We will consider the sensitivity of a LS neutrino detector with a spherical fiducial volume 
Vo ~ 1000 m^ and 2135 live-days of data (these are roughly the specifications of KamLAND). 
We estimate the neutrino detector's sensitivity utilizing the procedure outhned in section 
II. The density of the liquid scintillator is taken to be 80% that of water. Following [31] . 
we define as "fully-contained" an electron/positron event starting within the detector with 
at least 10 radiation lengths (~ 4.3 m) contained within the detector. Furthermore, the 
lepton event must point back to the sun within a half-angle 6cone = 20°-^10 GeV /Ey, and 
the energy of the neutrino must be obey E^ > 1.5 GeV. We then find [21] 

f 2 1 

/ dVr]{r) X e{r,z) ~ x - x -Vq, (21) 
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where the factor 2/3 is the fraction of lepton events which will be within Ocone-, and the 
factor 1/2 is the fraction of the fiducial volume which will yield fully-contained events. It 
was shown in [31j that, using the background estimate of [59], one would expect less than 
5 electron/positron events satisfying these cuts arising from atmospheric neutrinos during 
the specified runtime. We will thus consider a model which would produce 10 signal events 
arising from dark matter annihilating in the sun as being excludable. 

To determine the electron (anti-) neutrino flux at Kamioka, we have used the Wimp Event 
routine, run on the Hawaii Open Supercomputing Center cluster. Of the 10^ dark matter 
annihilations which were simulated (as described in section 111) for each annihilation channel 
and value of mx, 2 x 10^ were used to compute the neutrino spectra at the detector. The 
effect of neutrino propagation through the earth typically suppresses {Nz) by ~ 25 — 50% 
(depending on the mass and the annihilation channel). 

For CDMS, we will roughly estimate the sensitivity to dark matter with long-range inter- 
actions from their published bounds [7] on dark matter with elastic contact isospin-invariant 
interactions. The bound CDMS can place on C can be related to its bounds on cigj by the 
relation 

<^ ,^ ml r du[f(u)/u]w^E-lf^-''dER\FGeiER)\^ 

f^bound _ p(bound) '"-Ge Ju^in LJ \ J / J max J Ethr ^' ' /ooN 



IZ. [fiu)/u]w^E-l^ dEn \Fa.{EnW/E],' 



where u is the velocity of a dark matter particle far from the sun, and w is the is the 
velocity of the same particle once it has reached the surface of the earth. When it reaches 
the surface of the earth, the kinetic energy of the particle has increased by an amount 
equal to the change in the gravitational potential energy. The change in the gravitational 
potential energy due to the sun and the earth (Kun and Vearth, respectively) can be written 
as AVsun,earth = -(l/2)mxf Ln,eari/i5 whcre Vgun ~ 42.1 km/s is the escape velocity of the 
sun at the radius of the earth's orbit, and Vearth ~ 11-2 km/s is the escape velocity of the 
earth at the surface of the earth. Using the relation AEkmetic = —AVsun — AVearth, we find 
w = (u2 + i;2^„ + t;2^^jl. 

The maximum recoil energy which can be transferred to a germanium nucleus is Emax = 
2m\mGew'^ / {.'f^x + ^GeY- We assume a threshold energy Ethr = 2 keV |7]. Umin is the 
minimum dark matter velocity (far from the sun) such that scattering with Eji > Ethr is 
kinematically possible, and is given by the expression u^^^ = max[(mGe-E't/ir/2m^) — — 

'^earthy 0] • 

We again assume a Gaussian form factor Fa{Epi)] for the recoil energy range of interest, 
the Gaussian form factor for germanium differs from the Helm form factor [60j by at most 
6%. We will assume a Maxwell-Boltzmann velocity distribution with v = 270 km/s, and 
that the galactic escape velocity is 600 km/s. We will also assume a constant efficiency for 
events with recoil energy greater than the threshold energy to appear in the CDMS low- 
energy analysis band. Due to this assumption, the result shown here should be regarded as 
only an estimate of the sensitivity CDMS could obtain with present data to dark matter 
models with long-range interactions. 

The estimated sensitivity of CDMS and a 1 kT liquid scintillation neutrino detector are 
plotted in figure |6] For the LS neutrino detector, we assume 2135 live-days of data, and 
assume that dark matter annihilates exclusively to either rf, bb, cc, gg or (with equal 
coupling to all three neutrino flavors). In [3T], it was shown that the sensitivity of CDMS to 
10 GeV dark matter with isospin-invariant elastic contact interactions is roughly an order 
of magnitude greater than that of a 1 kT LS detector. Our calculation of the relative 
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FIG. 6: Sensitivity to C of CDMS and a 1 kT LS detector (2135 live days of data) for low-mass 
dark matter with isospin-invariant elastic long-range interactions. LS detector sensitivity is shown 
assuming annihilation to either the r, h, c, g and v (flavor- independent) channels. 

sensitivities of CDMS and a 1 kT LS detector to dark matter with long-range interactions 
bears out our original estimate of a roughly 10^ — 10^ relative enhancement in sensitivity 
for the LS detector. Note that, for the models to which KamLAND would be sensitive, 
the sun would be in equilibrium (see figure |5]) even if the annihilation cross-section were 
significantly smaller than 1 pb (assuming standard astrophysical assumptions). If the sun 
is not in equilibrium as a result of deviations from these assumptions, then the constraints 
which would be possible from neutrino detectors would be significantly suppressed. 

We will not attempt a quantitative estimate of the sensitivity of XENONIOO to dark 
matter with long-range interactions. XENONlOO's recoil energy threshold is defined in terms 
of scintillation photoelectrons; the detector's scintillation response to recoil energy (£e//) 
is not measured for low recoil energies. Moreover, bounds from XENONIOO are generated 
assuming that the number of photoelectrons is determined by a Poisson distribution. Even 
some low energy recoils can thus produce enough scintillation photoelectrons to exceed 
the threshold. As a result of the uncertainties in the detector response at low recoil, an 
attempt to estimate the event rate expected at XENONIOO for dark matter with long-range 
interactions is beyond the scope of this work. We will simply note that the recoil energy 
range for which XENONIOO is sensitive is at best comparable to that of CDMS, while 
a xenon nucleus is roughly twice as heavy as that of germanium. This implies that the 
sensitivity of XENONIOO relative to CDMS will be suppressed by roughly a factor of 4 for 
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the case of long-range interactions. 

Recent hints of low-mass dark matter have potentially been seen by the DAMA [1], 
CoGeNT [2] and CRESST [3J experiments. Dark matter models with long-range interactions 
have been discussed as a possible way of reconciling the data from these experiments with 
the constraints from other direct detection experiments [SI1E2]- Long-range interactions can 
affect not only the magnitude of the overall excesses seen by CoGeNT and CRESST, but can 
also affect the modulation seen by DAMA and CoGeNT. Although models with long-range 
interactions can provide a better fit to the overall excesses, they sometimes provide a worse 
fit to the modulation signals. We will not attempt to define a region of parameter-space 
for dark matter with long-range interactions which could match the DAMA, CoGeNT or 
CRESST data. This would require a detailed matching of the expected event spectrum with 
that observed by the experiments, which is beyond the scope of this work (and perhaps 
premature, given the issues raised in [IS])- However, since CoGeNT also uses germanium as 
the target material, one would expect a 1 kT LS detector to be easily sensitive to dark matter 
models with long-range interactions which could potentially explain the data of CoGeNT. 
Moreover, one should note that the low-mass CRESST region is consistent with scattering 
from both oxygen and calcium. Given the difference in mass, one would expect long-range 
interactions with calcium to be suppressed by roughly a factor 4 relative to oxygen, as 
compared to the case of contact interactions. 

Finally, we can consider the sensitivity of LBNE (Long-Baseline Neutrino Experiment). 
We will assume the detector target material is liquid argon (configuration 2) [HS], with a 
total fiducial volume of roughly 51 kT. Liquid argon-based neutrino detectors are expected 
to have very good event reconstruction; we will assume that liquid argon detectors permit a 
reconstruction of charged lepton flavor, energy and direction with at least the same resolution 
as liquid scintillation detectors. We then find 



where the factor 2/3 again arises from the fraction of charged lepton events which would 
point back to the sun within angle 9 cone- For a detector as large as LBNE, almost the entire 
fiducial volume can produce fully-contained events. We then see that the sensitivity which 
KamLAND could obtain with its 2135 day data set could be obtained by LBNE with only 
~ 17 days of data. 

Note that the possibility of dominant annihilation to leptonic channels is not inconsis- 
tent with dark matter-nucleus scattering which is large enough to be probed by neutrino 
detectors. For example, it could be that dark matter-quark scattering is mediated by an ef- 
fective operator which permits velocity- independent, spin-independent scattering, but does 
not permit s— wave annihilation (an example of such an operator is XXqq). In this case, 
the dark matter-nucleus scattering cross-section could be reasonably large, while the the 
cross-section for dark matter to annihilate to quarks would be u^-suppressed. If dark matter 
coupled to leptons through an operator which permitted s— wave annihilation (an example 
of such an operator would be X'-y^X f'-y^ f , if the dark matter were a Dirac fermion), then 
the dark matter would mostly annihilate to leptons. 




(23) 
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VII. CONCLUSIONS 



We have computed the capture rates and neutrino spectra which are relevant for neutrino- 
based searches for low-mass dark matter in the sun. The neutrino spectra are presented at 
a distance of 1 AU from the sun, accounting for matter effects in the sun, and vacuum 
oscillations (assuming a normal hierarchy and 613 = 10°). The capture rates have been 
found assuming either elastic contact, elastic long-range, or inelastic contact interactions. 
These are the tools required for a neutrino detector to search for dark matter annihilating 
in the sun. 

As an application of these tools, we plot the sensitivity of a 1 kT LS detector, with 2135 
days of data, to low-mass dark matter with isospin-invariant elastic long-range interactions 
with Standard Model nucleons. We have found that neutrino detectors have a greatly 
enhanced sensitivity to dark matter with long-range interactions, relative to leading direct 
detection experiments such as CDMS. This enhancement is readily understood; in the case 
of long-range interactions, the scattering matrix element is inversely proportional to = 
2mAER. Scattering rates in detectors with heavy targets, such as germanium and xenon, 
are heavily suppressed. Dark matter capture in the sun involves scattering from low-mass 
targets such as hydrogen and helium, implying that these scattering rates will see a relative 
enhancement. A LS neutrino detector with the exposure already available to KamLAND 
could have a sensitivity up to 2 orders of magnitude greater than that of CDMS. LBNE 
(with a 51 kT liquid argon target) could achieve similar sensitivity with roughly 17 days of 
data. 

We have also found that low-mass dark matter with inelastic contact interactions can be 
probed by neutrino detectors even for 5mx ~ 50 keV. This implies that neutrino detectors 
can be sensitive to inelastic dark matter models which are more difficult to probe on earth, 
because gravitational infall allows inelastic scattering in the sun for models where inelastic 
scattering would not be kinematically possible on earth. 

The choice ^13 = 10° is consistent with recent data from the Daya Bay experiment [53] . 
The neutrino spectrum is slightly different from the ^13 = 0° case, with the difference 
most noticeable in the case of annihilation entirely to neutrinos. For searches involving 
upward-going leptons, there will also be a modification to the neutrino spectrum due to 
passage through the earth. This effect will depend on the location of the detector; for 
any particular detector, one can obtain the appropriate neutrino spectra by running the 
WimpEvent program, inputting the data files for the neutrino spectrum at 1 AU found at 
http : / / w ww . phys .hawaii . edu/~superk/post/ spectrum, 

It is worth noting that a direct detection experiment with a target molecule containing 
hydrogen would also be expected to have enhanced sensitivity to dark matter with long- 
range interactions. Gaseous time projection chambers (such as DRIFT [64j, DMTPC [03] , 
MIMAC [67j and NEWAGE [68j) using hydrocarbon targets may be well-suited for 
this type of search. 

Specific dark matter models with long-range interactions may have solar capture rates 
that are enhanced by collective effects, such as multiple scattering. Neutrino searches thus 
have enhanced sensitivity to such models, and current data may already provide tight con- 
straints. It would be interesting to consider such models in more detail. 



17 



VIII. ACKNOWLEDGMENTS 



We gratefully acknowledge K. Choi, D. Marfatia, M. Sakai, P. Subramoney and S. Vahsen 
for useful discussions. We also thank the Hawaii Open Supercomputing Center. This work 
is supported in part by the Department of Energy under Grant DE-FG02-04ER41291. 

Appendix A: Neutrino Spectra 

In this appendix, we plot the neutrino and anti-neutrino spectra at a distance 1 AU from 
the sun, assuming dark matter annihilation exclusively to either the bb or rf channels. Each 
spectrum was generated by simulating 10^ annihilations using the method described in the 
text. Here, z = E^/mx, and neutrino oscillations are generated assuming ^13 = 10°. 
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FIG. 7: Neutrino spectra (left panels) and anti-neutrino spectra (right panels) at 1 AU for dark 
matter annihilation to the bb channel. The spectra for Ve{i>e)-, ^'^(j^/i), and ^'r(^'r) ai'e shown in red, 
green, and blue, respectively. Spectra are shown for mx = 6,8 GeV. 
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FIG. 8: Neutrino spectra (left panels) and anti-neutrino spectra (right panels) at 1 AU for dark 
matter annihilation to the rf channel. The spectra for f^(i/^), and z^r('^r) are shown in red, 

green, and blue, respectively. Spectra are shown for mx = 4, 6, 8 GeV. 
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